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Abstract. We present the gravitational wave signatures for a suite of axisynimetric 
fSJ core collapse supernova models with progenitors masses between 12 and 25 M©. 

^ These models are distinguished by the fact they explode and contain essential physics 

(in particular, multi-frequency neutrino transport and general relativity) needed for 
r — a more realistic description. Thus, we are able to compute complete waveforms 

^^ {i.e., through explosion) based on non-parameterized, first-principles models. This 

l/~) is essential if the waveform amplitudes and time scales are to be computed more 

^^ precisely. Fourier decomposition shows that the gravitational wave signals we predict 

should be observable by AdvLIGO across the range of progenitors considered here. The 
fundamental limitation of these models is in their imposition of axisymmetry. Further 
progress will require counterpart three-dimensional models. 
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1. Introduction 

Core collapse supernovae are among the sources that will produce gravitational 
waves (GWs) detectable by GW observatories around the globe. In particular, a 
Galactic supernova event will likely produce a signal well within Advanced LIGO's [1] 
bandpass across a broad range of frequencies. Gravitational waves from supernovae 
will arise from a variety of phenomena given their multidimensional, multi-physics 
character. These phenomena include fluid instabilities in the proto-neutron star, 
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neutrino-driven convection beneath the supernova shock wave, the standing accretion 
shock instabihty (SASI), deceleration at an aspherical shock, and aspherical neutrino 
emission. Obviously, two-, and ultimately three-, dimensional models will be required 
to capture the GW emission from such phenomena. Moreover, the explosion dynamics, 
and ultimately the computation of the GW emissions as well, will require sufficient 
realism in the treatment of core collapse supernova multi-physics. For a comprehensive 
survey of the field see the recent review by Ott [2] and references therein. 

2. The CHIMERA Code 

Our CHIMERA code has five major components: hydrodynamics, neutrino transport, 
self-gravity, a nuclear equation of state, and a nuclear reaction network (see [3] 
for details). The hydrodynamics is evolved via a Godunov finite- volume scheme — 
specifically, a Lagrangian remap implementation of the Piecewise Parabolic Method 
(PPM) [4]. Neutrino transport along our radial rays is computed by means of multigroup 
fiux-limited diffusion (in the "ray-by-ray-plus" approximation [5]), with a fiux limiter 
that has been tuned to reproduce Boltzmann transport results to within a few percent 
[6]. A spectral Poisson solver is used to determine the gravitational field [7], with general 
relativistic (GR) corrections to the spherical component [8]. Details of the integration of 
this gravitational framework with the hydrodynamics will be given in [9] . The equation 
of state (EOS) of Lattimer & Swesty [10] is currently employed for matter in NSE above 
1.7 X 10^ g/cm . Below this density, matter in NSE is described by 4 species (neutrons, 
protons, helium, and a representative heavy nucleus) in a corrected and improved version 
of the Cooperstein EOS [11], extended to regions where the composition is determined 
externally by a reaction network. Continuity with the LS-EoS is achieved by use of a 
common electron-positrion EOS (a revised and extended version of that in [11]) and 
establishment of a common zero for the mass energy. For regions not in NSE, an EOS 
with a nuclear component consisting of 14 a-particle nuclei from ^He to ^°Zn, protons, 
neutrons, and an iron-like nucleus is used. An electron-positron EOS with arbitrary 
degeneracy and degree of relativity spans the entire density-temperature regime of 
interest. The nuclear composition in the non-NSE regions of these models is evolved by 
the thermonuclear reaction network of Hix & Thielemann [12]. While Eulerian schemes 
are preferred for regions with violent turbulence, they have a disadvantage: the history 
of field variables for a given parcel of material, crucial for nucleosynthesis, is lost. To 
compensate for this loss, and to allow post-processed nuclear network computations, the 
tracer (or test) particle method [13] has been implemented in CHIMERA. The tracer 
particles are equally distributed on the spherical grid (40 particles/row x 125 rows) at the 
pre-coUapse phase and follow the fiow in the course of the Eulerian simulation, recording 
their temperature and density history by interpolating the corresponding quantities from 
the underlying Eulerian grid [13]. Each particle is assigned a constant mass (1/5000 of 
the progenitor mass) and the GW signal it produces is calculated taking the quadrupole 
integral. Comparing the GW corresponding to a given group of tracers with the signal 
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Figure 1. Left: Entropy distribution at 244 ms after bounce for the 15 M0 model. 
A large, low-entropy (blue-green) accretion funnel at an angle quasi-orthognal to the 
symmetry axis and high-entropy (yellow-orange-red) outflows below the shock, along 
the symmetry axis, are evident. Right: Shock radius as a function of time for three 
regions: the north pole (solid blue), the equatorial plane (dotted black), and the south 
pole (dashed red). 



produced by the bulk matter motion allows us to identify what part of the fluid generates 
a specific GW feature. 

3. Gravitational Wave Extraction 



Stochastic matter motion and anisotropic neutrino emission during the explosion 
generate GWs. The transverse-tracefree part (TT) of the gravitational strain is written 
as 

1 m=2 / r \ 2 
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where the mass quadrupole (as a function of retarded time) is computed by 
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with Too the corresponding component of the linearized stress-energy tensor, and 
P^{9, 0) the spherical harmonics. In the weak-field case, we approximate Too — p, where 
p is the rest-mass density. Following the Finn-Evans approach [14] to reduce the second 
time derivative A2ra = ^hm = (i^2m and using the continuity equation (Blanchet et 
al [15]), we can calculate A^2m as in Equation (34) in [14]. In axisymmetric cases, N20 
is the only non-null component, and we evaluate its time derivative numerically. The 
wave amplitude is related to the dimensionless gravitational strain, h+, by 
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where r is the distance to the source and 6 is the angle between the symmetry axis and 
the observer's line of sight (we will assume sin^ ^ = 1). To compute the GWs produced 
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by anisotropic axisymmetric neutrino emission, we use Epstein's [16] and Miiller & 
Janka's formalism [17]: 

4^ r... r.n>..n>.dL{e',t') 



hr = ^ f dt' r de'^{e' 

c^r Jo Jo 



dn' ' 

with \l/(^) given in [18]. The direction-dependent differential neutrino luminosity, 
dL/dQ, is calculated at the outermost radial grid zone. In order to determine the 
detectabihty of the GWs, we calculate the characteristic GW strain for a given frequency 
/ [19] using 



1 / 2 G dEcw if) , dEawif) c' i^Tiff ~ 2 

where dEcw (/) /df is the GW energy spectrum and A20 (/) is the Fourier transform of 

A20{t). 

4. Gravitational Waveforms 

We performed axisymmetric two-dimensional simulations beginning with 12, 15, and 
25 Mq non-rotating progenitors [20] and resolutions of 256 (adaptive) radial and 256 
angular zones. The radial grid ranges from km to 1.88-10^ km. The left panel of Figure 
1 shows an entropy snapshot of the 15 Mq model. Successful explosions are obtained in 
all cases, with the longest running model (the 25 M© model) having an explosion energy 
of 0.7 B (and still growing) 1.2 seconds after bounce. Details are provided in [21, 3]. 

While the GW emissions we predict differ in detail from model to model, a clear 
GW signature, composed of four parts (left column of Figure 2) emerges: 1) A prompt 
signal: an initial and relatively weak signal that starts at bounce and ends at between 
50 and 75 ms post-bounce. 2) A quiescent stage that immediately follows the prompt 
signal and ends somewhere between 125 ms and 175 ms after bounce. 3) A strong 
signal, which follows the quiescent stage and is the most energetic part of the GW 
signal. This stage ends somewhere between 350 ms and 450 ms after bounce. 4) A 
tail, which starts before the end of the strong signal at about 300 ms after bounce and 
consists of a slow increase in rh^. This tail continues to rise at the end of our runs. 

Waveforms covering the first three of four phases (prior to explosion) have been 
computed by Marek et al [22], and waveforms covering all four phases and based on 
parameterized explosions were reported in the work of Murphy et al [23]. The overall 
qualitative character of the GW signatures shown in [23] reflect what is shown in Figure 
2. The work presented here takes the natural next step beyond this earlier, foundational 
work. A more precise prediction of the GW amplitudes and the timescales associated 
with each of the four phases requires a non-parameterized approach. Even in the case 
of a non-parameterized approach, prior to evidence of an explosion it is difficult to 
assess whether or not the amplitudes and timescales are well determined. Thus, the 
non-parameterized explosion models studied here enable us to predict all four phases of 
the GW emission and their amplitudes and timescales with some confidence. 
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Figure 2. The left column shows the GW strain times the distance to the observer 
vs. post-bounce time for non-rotating progenitors of 12, 15, and 25 Mq. The signal is 
split into matter- (red-solid) and neutrino-generated (blue-dashed) signals. Note that 
the scales are different for these two signals. The insets show the first 70 ms after 
bounce. The right column shows the corresponding characteristic strain for both the 
matter (red) and the total (black) signals, compared to the AdvLIGO sensitivity curve. 
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Figure 3. Contributions to the matter-generated GW signal from two different regions 
for the 15 Mq model: the PNS (r < 30 km) and the region above the PNS (r > 30km). 
The latter includes the region of neutrino-driven convection, the SASI, and the shock. 



The prompt signal is generated by two independent phenomena: prompt convection 
inside the PNS generates a high-frequency signal that is superimposed on a lower- 
frequency component, seen in the insets of Figure 2. There is a hint of this in the inset 
of Figure 3, where the signal for our 15 Mq run has been split into the contributions 
from two different regions, but it is in the tracer analysis of Figure 4 (right) where this 
becomes evident. The matter-generated GW (solid red) is closely tracked by the GW 
generated by the infalling tracer particles deflected by the shock (dashed blue), some 
of which are shown in the left panel. The low-frequency signal from 20 ms to 60 ms 
after bounce originates at the shock radius, which is at ~100 km at this time and well 
outside the PNS. In the past, authors attributed the prompt signal to convection only 
[22, 23]. 

The quiescent stage corresponds to the period after prompt convection has ceased 
and prior to the development of neutrino-driven convection and the SASI. It is followed 
by a strong signal produced by the development of both. The strong signal is dominated 
by SASI-induced funnels impinging on the PNS surface. It shows evidence of two 
components (also described in [22, 23]): The low-frequency component arises from the 
modulations in the shock radius as the SASI develops and evolves. The right panel of 
Fig. 1 shows the first cycle of this modulation at 1 75ms with the maximum (minimum) 
of the north (south) pole radius and at about 210ms with the reverse situation. The 
high-frequency component is generated when the SASI-induced accretion flows strike 
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Figure 4. Le/i; Trajectories of the tracer particles. It is sliown tlie clear deflections 
of infalling particles through the shock that collectively produce low-frequency high- 
amplitude component of the GW signal shown on the right panel. Right: Comparison 
between the matter signal (solid red) and signal calculated using the tracers (dashed 
blue). Both panels correspond to our 15 Mq simulation. 



the PNS (Fig. 1 left). The shock modulations affect the kinetic energy of the accretion 
flows and, consequently, the amplitude of the GWs generated when these flows hit the 
PNS. Hence the high-frequency modulations are beneath a low-frequency envelope. 

All of our GW signals end with a slowly increasing tail, which reflects the 
gravitational memory associated with accelerations at the prolate outgoing shock (see 
also Fig. 5 in [23]). The tail continues to rise at the end of our runs because the 
explosions are still developing and strengthening. The 15 M0 model GW is shown in 
Figure 3, where the explosion starts at ~ 300 — 350 ms after bounce, and by 400 ms the 
signal from the PNS has largely ceased. 

Focusing now on /ichar (Figure 2 right), it is important to note that the peak at 
~ 700 — 800 Hz is associated with the high-frequency component of r/i+, which in 
turn is associated with the downflows hitting the PNS surface, as discussed above. 
A precise association of the signal at lower frequencies with phenomena in the post- 
bounce dynamics will require a detailed analysis using tracer particles and will be left 
to a subsequent paper [9]. The lower-frequency modulations (the envelope) in r/i+, 
which in turn are associated with the SASI-induced shock modulations will certainly 
be an important component of this lower frequency signal. Finally, we note that the 
possible dependence of the AdvLIGO-observable signal between 100 Hz and 700 Hz on 
the progenitor mass also requires further investigation and a detailed discussion, which 
will be presented in [9]. 

The amplitudes of the GWs from neutrino emission are negative from bounce 
to ~ 180 — 220 ms after bounce and then increase dramatically, becoming positive 
throughout the end of the simulation. The positive sign is consistent with a relative 
dominance of neutrino emission along the polar over the equatorial regions [18]. The 
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Figure 5. Energy emitted by GWs during the first 500 ms after bounce for all three 
models presented here. 



change in sign from negative to positive correlates with the formation of the funnel- 
hke downflows of dense matter, which increase neutrino opacities in the equatorial 
plane (orthogonal to the symmetry axis; see Figure 1). Note that the amplitude 
of the neutrino-generated GW signal is much larger than the matter-generated GW 
signal. However, these GWs have relatively low frequencies, and their contribution to 
the total characteristic strain is only significant at frequencies below 20 Hz (see also 
[17, 18, 24, 25]). 

Our GW predictions for the 15 M© case can be compared to those of Marek et al 
[22] given that both groups implement similar treatments of the neutrino transport and 
GR corrections to the gravitational field, and include essentially the same overall multi- 
physics in their models. The two groups are in agreement with regard to the time scales 
of the different (pre-explosion) GW phases, the amplitude of the prompt signal, and the 
peak in the GW spectrum at ~ 700 — 800 Hz. They differ, however, in their predictions 
for the amplitude of the strong signal, where our results are about twice as large. These 
differences will be investigated. They likely arise in part due to the different progenitors 
used, which in turn alters the time scale to explosion and, consequently, the amplitude 
of the GW signal in the strong-signal phase at any instant of time. 

A comparison with Murphy et al [23] is more difficult given that their models 
are parameterized (i.e., their neutrino luminosities are kept isotropic and constant 
throughout their simulations). As a result, we can expect the GW amplitudes, time 
scales, and frequency peaks to be different. Murphy et al [23] present a simple model 
for the GW strain during the strong-signal phase, where the GW amplitude, caused by 
the downdrafts hitting the PNS surface, is proportional to the downdraft frequency fp 
and velocity fp, which, in turn, depend on the compactness of the PNS. In light of this 
model, the larger amplitudes and characteristic frequencies we see might be attributed 
in part to our use of an effective GR potential and the soft Lattimer & Swesty EOS 
rather than a Newtonian potential and the stiff Shen EOS [26]. 
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The total emitted GW energy is shown in Figure 5. For the more massive 
progenitors, all of the GW energy is emitted between 200 ms and 400 ms after bounce. 
For the 12 M© case, the GW energy is emitted more slowly, consistent with the fact 
that the explosion in this case unfolds more slowly [3]. Our predictions are 20 to 50 
times larger than those of Murphy et al. [23] , but this is consistent with our waveforms 
having two to three times the amplitude and a higher frequency than the signals they 
predict. 

5. Summary and Conclusions 

We present gravitational waveforms computed in the context of 2D core collapse 
supernova simulations performed with the CHIMERA code for non-rotating 12, 15, 
and 25 Mq progenitors. We calculate the contribution to the signals produced by both 
baryonic matter motion and anisotropic neutrino emission up to 530 ms after bounce 
for all three progenitors. Given the development of non-parameterized explosions in 
our models, we are able to compute the waveforms through explosion and to determine 
more precisely the pre-explosion amplitudes and timescales. Given our use of tracer 
particles, we are able to decompose the GW signatures and determine which phenomena 
contribute to specific components of the waveforms. This allowed us to identify 
an additional source for the prompt signal (in the past solely attributed to prompt 
convection) the deflection of infalling matter through the shock. 

Our waveforms exhibit a characteristic signature. Namely, the signal develops in 
four stages. There is 1) a relatively short and weak prompt signal, 2) a quiescent 
stage, 3) a strong signal where most of the GW energy is emitted and, 4) a slowly 
increasing tail. We predict signatures with sufficient strength to be readily observable 
by Advanced LIGO for a Galactic event, and the peak in the observable spectrum stems 
from the accretion downfiows driven by the SASI. 

The results presented here are preliminary: a new set of 2D simulations performed 
with an enhanced version of our CHIMERA code is currently ongoing [9]. However, 
while a number of approximations are made in the CHIMERA code, the more important 
limitation in these models is their restriction to axisymmetry. Three-dimensional models 
are required. We anticipate that the greatest change to our gravitational waveform 
predictions in moving to 3D will be in the phase-4 tail. Prolate explosions are often seen 
in axisymmetric simulations, where artificial boundary conditions must be imposed that 
prevent the turnover of material along the symmetry axis. With axisymmetry removed, 
we expect a significant change in the evolution of the explosion tail: in its magnitude, 
and perhaps even in its sign. And no doubt there will be quantitative changes to the 
amplitudes and timescales associated with earlier phases, particularly in the strong 
signal arising from the SASI motions. In 3D, the SASI will likely be dominated by 
spiraling fiows [27], fundamentally different than the sloshing modes that dominate in 
the axisymmetric case. This will in turn alter the waveforms in the final pre-explosion 
phase. This has already been demonstrated by Kotake et al [25] in 3D parameterized 
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studies. 3D simulations with all of the physics documented here are ongoing, and we 
look forward to reporting on their GW signatures in the near future. 
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